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SOQ USER GUIDE UPDATES 


June 1980 Updates to SOQ80128 


INTRODUCTION 


This document defines the changes made to the SOQ code (S0Q80128) 
between January and June of 1980. The changes either correct short- 
comings found in the code or, more usually, document the increased 
capability being continually built into the code. The SOQ code is 
raintained as S0Q80128 June PL,ID = ARLOJRA as a NOS/BE-1 CDC update 
format file. 


UPDATES 

1. ¥*ID FIXZRN 

This update redefines the coeffiencets to be input to the Zernike 
subroutine. This new convention is more physically meaningful in that, 
at least for lower orders, the coefficients are in waves. For example, 


to impose one wave peak to peak of defocus (P,) on a beam, one would 
input P(4)=1. The phase applied is now: 


o(1,J) = EP mz (1,J) 


The subroutine affected is ZERN. This update does not effect the rest 
of the code. 


2. ¥*ID FIXJTR 

This update ensures a correct definition of DF in subroutine JITRBG 
since when JITRBG is called from subroutine QUAL, the X-coordinate 
array contains RA/D coordinates, not the spatial coordinates. 

Only one line of the code is affected by this update. 


3. *ID ROTZRN 


Due to different coordinate system orientations for data, it became 
necessary to allow for this variation within subroutine ZERN. 


Define the data x and y coordinates to be XROT and YROT, and the SOQ 
x and y coordinates to be XIN and YIN. The rotation angle is then 
defined to be © (in radians). 


<a 
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COSROT = COS(@) 
SINROT = SIN(0) 


XROT XIN x COSROT 


+ YIN x SINROT 


YROT ~XIN x SINROT + YIN x COSROT 


Application of Zernike polynomials to and SOQ point located at 


(XIN, YIN) would then be 
bility of axis flips are 
FLIPX or FLIPY not equal 
include FLIPX, FLIPY and 


calculated using Z(XROT, YROT). The possi- 
also accounted for and are flagged by 

to zero. Namelist ZERNS is modified to 

the rotation angle (in degrees) ZTHETA. No 


common was modified. This update modified only subroutines GDL and 


ZERN. 
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SECTION V APPENDICES 
APPENDIX A 
SOQ USERS GUIDE UPDATES 
JANUARY 1977 TO JUNE 1978 


INTRODUCTION 


This appendix documents those changes made to the initial SOQ code between 
January 1977 and June 1978. The changes incorporated in the code are those that 
have become generally useful for the physical optics simulation problems which 
have been solved using the SOQ code. The Users Guide Updates are also prepared 
to clarify and correct the initial description of the SOQ code, as documented 
and delivered to AFWL on 1 March 1978, in the Preliminary SOQ Users Guide. 

This document supersedes previous written material on SOQ code documentation. 


The organization of the SOQ Users Guide Updates is 


SECTION AI New Subroutines 
1. Theory 
2. FORTRAN Updates 
SECTION AIT Code Changes/Corrections 
1. Theory/Reason for Correction 
2. FORTRAN Updates 
SECTION AIITI Users Guide Corrections 


SECTION AIV SOQ Code Access 


AI. NEW SUBROUTINES 
1. THEORY 


a. Beam jitter -- Relative motion between optical elements, such as 
mirrors, induces time varying positional displacement of the optical field. 
The typical term for this phenomenon is beam jitter, and the principle effect 
is to broaden the sime-averaged effective beam illumination area, while 


reducing the time averaged intensity. 


Beam jitter is both a near field and far field concern. Jitter in optical 


trains can overload apertures or cause energy deposition on areas outside the 


normal propagation path as well as cause a deterioration in the peak on-axis 


irradiance and integrated spot power. 
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The time-averaged effect of beam jitter may be modeled as the convolution 


of the intensity profile with an appropriate probability density function (PDF) 


for the jitter statistics. The current SOQ model assumes that the jitter PDF 


is Gaussian with known mean and variance. The model allows the user to speciry 


the Gaussian parameters and, for the selected beam jitter analysis location, to 


determine the near field and/or far field effect of beam jitter. 


The following is a brief description of the analytical and SOQ Fortran 


implementation of beam jitter calculations: 


b, Relevant formalism -- The effect on the beam may be found by convolu- 


tion of the Gaussian jitter probability density function with the SOQ predicted 


intensity distribution: 


J 
1’ (x,y) - ff: Ch y)J(x-x’, y-y! ) dx’ dy’ (Al) 
-00 


The 1-D Fourier transform of the Gaussian function looks like: 


2 pee rRe 


H(f) 


0 


FORTRAN UPDATES 


Each assumes that the jitter 


The jitter model can be called in two ways. 
variance is the product of a jitter angle and the propagation distance from the 


jitter source. 
(A2) 


= Jitter angle (lo, in microradians) 


= Distance from jitter source (in cm) 


When the far field model is called from QUAL, the jitter angle has been 
incorporated into namelist QLOT while the propagation distance is the focal 
length found in QUAL. The jittered intensity is returned to array CU as a 
phaseless field so it can be plotted, or written to a permanent file. 


The other method of activating the jitter model is to call the near field 
jitter model from GDL with IFLOW = 23. For this model both angle and jitter 
distance are entered in namelist JITTER. 

Namelists modified: 
Far Field 
QLOT: SIGANG (rad) is added to specify the jitter angle 


Near Field 
JITTER: Contains - 
JITANG (urad) 
JITDIS - Jitter distance 


ATI. CODE CHANGES/CORRECTIONS 


THEORY/REASON FOR CORRECTIONS 


a. Bare resonator calculations -- The SOQ resonator/optical-train cal- 
culation code may be used to simulate, in Cartesian coordinates, bare reso- 
nators. This added option is frequently used in the initial simulation studies 
of a resonator or a class of resonators. 


The bare resonator optical configuration may be compared to its geometric 


counterpart using the SOQ code by simply invoking the IBARE option and asso- 


ciated updates now contained in the fundamental code. The fundamental approach 


in bare resonator calculations on the SOQ code is to allow the user to use the 
same input and code for bare, semibare and loaded cavity calculations. Various 
options under the bare cavity calculations have been incorporated and are now 
described as input values for IBARE in Namelist START. 
IBARE = 0 (Default) 
Loaded cavity calculations are performed as usual following the stan- 
dard input which the user has supplied. 
IBARE = 1 
Using the same input, the user will now perform bare cavity calcula- 
tions in which the resonator is normalized to 1W of circulating power. 
Mirrors are defaulted to have 100 percent reflectivity, and no power 
dependent or flux dependent distortion. The SOQ output is modified 
“to printout the resonator eigenvalue. 


IBARE = 2 
Semibare resonator calculations are performed in which the user can 
perform bare resonator calculations that include optical aberrations 
generated by a flowing saturable gain medium. These aberrations may 
strongly effect mode shape/phase. This option provides a convenient 
method of studying their perturbational effect on the bare cavity mode. 


For the semibare option, an additional update has been included in which 
the namelist MIRROR user may specify the desired power at each of the resonator 
mirrors. This allows the user of the semibare updates to apply mirror distor- 
tions as though the bare cavity mode had a significant power level. Specifica- 
tion of the DESIPW value in namelist MIRROR to some value other than 0.0 will 
cause the field incident on the mirror to be scaled to that power specified by 


the numerical value of DESIPW. Appropriate mirror distortions will be applied 
at the desired power level. The field leaving the mirror will then be rescaled 
to its incident power level. Subsequent calculations are done as specificed by 


the user typical namelist input. 


An additional variation is allowed in which the parameter FLAG of cavity 
input namelist CAVTY2 can be used to execute a resonator with loaded gain, but 
no fixed phase perturbation in cavity. The input would correspond to FLAG = 0; 
IBARE = 0. Usual loaded resonator calculations are performed with mirror dis- 
tortions as specified by the user. 


All of the above variations of cavity/resonator calculations may be run 


from the standard loaded cavity input. 


At the rear of this section are Fortran listings of the code updates which 
have been included in the basic Cycle III SOQ code previously documented. 


b. Mirror non-normal incidence angle -- In many optical train calcula- 
tions the propagating optical field is incident on the mirrors in a nonnormal 
manner. Since, in general, the mirror surface may have a spherical figure, the 


field leaving the mirror will exhibit phase front aberrations introduced by non- 


normal incidence of the field on the curved surface. 


The SOQ MIRROR subroutine has been modified to incorporate an astigmatic 
aberration due to the nonnormal incidence on a spherical surface. The follow- 
ing is a brief description of the generation of the astigmatic aberration 
applied. 


Astigmatism in Resonator: 


General astigmatism is introduced when a wavefront is incident on a spheri- 
cal (parabolic) surface in a nonnormal manner. This aberration occurs at each 


spherically-distorted turning flat, for example. 


z 
y x - Sagittal Plane 
y - Tangential Plane @) 
Tg ts! ROsvG 
Ss 3”, : Ro 4 = Incident angle 
R, S #2 Object distance 
imi Cs Se 
Limite §% 2 cos¢ R_ = Mirror surface curvature 
. Cc os 
(spherical) 
Ae 2 ot x? S, = Sagittal plane effective 
SOQ, r 285, curvature 


Thus 28° is the resultant phase curvature being imposed on the beam. A 
cylindrical mirror can be used to model this with 


; R. = 2s° z -V2 R. (neg since Ro is convex) for 9 = 45 
he s/ 
(A3) 


Therefore, to represent the astigmatism introduced in the x-plane by a 
spherically-distorted turning flat, a cylindrical mirror is employed with a 


radius of curvature 


Sn -V2 R, 
Ss 
Re is the power induced radius of 
curvature which is input or deter- 
: mined by the SOQ code. 
© Similarly, the tangential nlane is described by 


2 
us y 
*8s0q, 7 x (42) 
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-R,cos¢ 
es, ean (A4) 
= —& for o= 45° 
2v% 
The new mirror subroutine including astigmatic effects has the form 
2r x? : 
8509 7 | 3er * Sr 
s T. 
R R_ cos 
Ss’ 2 —_ oc _ Si. = c 
s 2 cos ¢ rT 2 (AS) 


The only additional input change is to the MIRROR routine namelist which 
is expanded to include the variable PHIAST, the beam incidence angle in degrees 
(default is PHIAST 2 0). 


c. Beam rotation -- The mirror model has been updated to describe beam 
rotation introduced by optical elements. which are oriented in a skewed fashion. 
Many examples of this type of orientation are encountered in resonators and 
optical trains. The principle effect of skewed, or out-of-plane, orientation 
is to convolve or smooth the mirror distortion-induced aberrations over the 
total number of optical elements. 


Rotation of the beam is accomplished by analytically rotating the mirror 
with respect to the beam, rather than rotating the beam within the mesh and 
then applying the mirror. By rotating the mirror with respect to the beam 
two modeling advantages result: First, analytical rotation of the mirror with 
respect to the beam is accomplished with no interpolation loss of information. 
Second, since the rotation is analytical, computer time is saved by not having 


to evaluate the field numerically. The following describes the rotation 


equations used in the code. The following sketch shows a base and rotated system. 


y 


a = Rotation Angle 
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xX = X cosa + y sina (A6) 


y = -x sina + y cosa (A7) 
Then, 
2 
2 
q Here, 


(x,y) are the SOQ coordinates 


(x,y) are the transformed (rotated) coordinates 


The SOQ field is modified as 


; A9 
cu A exp (49) CULy (a9) 


OUT 


where A represents the completed transmittance effects included in mirror. 


The variable added to the SOQ MIRROR namelist input is PHIROT, which is the 
beam rotation angle in degrees. The default value is PHIROT = 0.0. 


2. FORTRAN UPDATES 


The attached printouts contain a listing of the updates, denoted as ROT, 
used to effect these changes. 


Litheal WOT 


@LOFNT wat 
“DELETE ClOaste.) 
ATOP (Dee) oXSCA 6G) eAHC (1402909) o FITLEI (20) oXUPANN (4) 6 
*ORLETF CLOFL A. 
DIMPNSTUN IN[9 (5626) oe Tal (99) oAHC (14020 09) eCFFL (16386) e LUSK (409) » 
®URLE TE CIOASTGL2 eCIOASTG.3 
DATA WANES eDOUTYOUINY oPH(AST ePHIHUr /500./ 
@PUVELETE CLOASTR. 
A DELTA SUISTR sIOUTY SULNY eH ANULS ePHIADTePMIROT NE SIP W 
®INSERT C1NaSTSG.S 
Cc MHTMOT TS Tree HEAM ROTATLI Ul ANOLE AT TraT STATI ANee DEG 
C NESTPwW TS TmME NESTRED POwEN LEVEL AT TeaT STATIUN 
®INGEWT C1NASTG YA 
AAC (1 %UeIMTRee) BON[ PUT 
. AAC (L@eIMIRe?) 2 NEST Pw 
2e OESTROw 2 AHC(Les{MIRoed} 
HELE TE ClOasNTGed 
; XK OTSTE eWANULS ePYOUT ORY INGMHTASTePMINT NEST Pw) 
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PINSERT MIRWOR, 2A 
PH TROTsPrIne 
*INSERT Cinasts.t2 
PHIRURS@PHIMUTOePI/ BU. 
PHIROT2Z0.0 
aR ITE (6086) PHI RU 
SINPRASIN(PRTROR) 
COSPRSsCUS (PMT RUR) 
*NELETE C10aSfG.15 
XPRMBX (J) OCOSPH OK (1) OSINER 
YPRMaek (J) OSINPR OK (1) *COUSPH 
PHASE SAK Y®( (XPRM2O2/HMSAG) © ( YHRMOC2/RMTAN) ) AK YONELL 
*VELETF C1HASTG.26 
APRMax (J) OCOSPHOX (1) OSINPR 
YPRMaek (5) OSINPWOk (1) PCUSPH 
PHASE SAKY® ( (KPRM®O2/HMSAG) © ( YPHM@OCs/OMTAN) } 
DELETE MIRROR .BS 
PHI Re (PHIASTOEPT)/180. 
RMSAGSRUC/COS(PHIR) 
RMTANSWOC*CDS (PHIR) 
PHI RORSePHIROTORI/140. 
PHlhOTsO.0 
SINPRSSIN(PHIRUR} 
COSPRSCOS (PHI RUN) 
*PNELETE MIPROR.DI 
XPRMBK (J) OCOSPREK (LT) OSI NER 
YPRMSek (J) OSINPROK (1) @COSPH 
PHIMIR SAR Y® ( (XPRM@O2/RMSAG) © ( YPM®O2/QMTAN) } 
CINSERT MIPROR.108 
WHITE (6046) PHIROR 
WRITE (46421) RMSAGeRMTAN 
-T-) FORMATI(2ZN40@ MIRROR RUTATIUN = %9612.50°¢R 00S) 


AIIT. USER'S GUIDE CORRECTIONS 
SUBROUTINE FUHS 


a. Purpose -- Subroutine FUHS is used to calculate the phase change 
due to heat release as the molecules in the lower laser level decay to the 
ground state. The FUHS modeling includes the assumption generally made for 
supersonic flow and assumes the heat release has only a small perturbative 
effect on the flow. 


b. Formulation -- The equations used here are based on those described 
by Biblarz and Fuhs, (Ref. 10) and-Fuhs, (Ref. 11). 


The usual continuity, momentum, and energy equations for steady flow 
with heat addition are used as the basis for the analysis: 


Continuity: v- (pu) = 0 


Do = 
Energy: Vepu h + & 


Momentum: 


= 


These are linearized assuming 


= 9 ’ a» 
Gei(+1u)+jyv 


Resulting in 


Continuity: 9, uw, + vy + Up’ = 0 (ALL) 


rey 
8 
x 


(Ys = uw etc) 


Momentum: 


Energy: 


The solution to these equations was first shown by Tsien and Bieloch, 


(Ref. 12), resulting in the following equations for a heat source q in super-~ 


sonic heat addition. 


ihe -(y - 1)q 6(x - 8y) 
2yp8 


(A15) 


», &-1) 6{x - By) (A16) 
as ay cca 
p= (r= 1Mg 8(x - By) | (y= a Sly) I(x) (A17) 


2a°B a“U 


Defines a Mach line 


Speed of Sound 


For a small volume, the heat addition is q = h{x,y) dx dy. The effects of all 


sources are then added; for example, 


Pees a - 
pt [ror 5(x-gy) dx dy 


s 


(y-1) ; 
=- “S54 f h(x=6y) sin yds (A18) 
rf 


where the integral is taken along a streamline (x = y) and siny = 1/M. 


s is related to x and y by 


xX = S COSH y = s sing (A19) 
By the Faltung theorem, for Fourier transforms, this can be written 
I’ (x,y)= P| F [tcx,¥) | °F [sen] (A20) 


The Fourier transform of the intensity is performed by the FFT, while the 


Fourier Transform for Gaussian density functions can be found analytically 


as 
22 
" me, exp Easy exp [ -2n0” (£7 + £)| (A21) 


From sampling theory, the discrete values of f. and fy can be found since 


A£ = + (A22) 


where 


L is the width of the SOQ calculation region (DCALC) 


£,) is then (I-1)4f. Recall from the discussion in FOURT, the DC value is 
returned in the first position and the last half of the transformed data are 
really negative frequency information shifted by one period, illustrated 


below in one dimension. 


Ax = Sampling rate in real space 


The equation for density change is, therefore, 
s 


Sk {s os { h(x,y) ee sin u ds 


2a°8 


(rs off dy’h(xyy’) 6(y = y’) T(x - x) 
aU 


The first term describes the compression waves along the streamlines due 
to heat addition, while the second describes the wake resulting from those 
compression waves. 


The heat release h(x,y) for a laser can be written: 


* ~ xX) 
h(x,y) = ¢ freee Rt 


NEP 


where t is the time constant for the depopulation of the lower laser level. 
If the depopulation were instantaneous (t + 0), the heat release would be 
proportional to the intensity, since every molecule emitting a photon would 
then immediately relax to the ground state with an accompanying increase in 
translational energy. It has been shown that the above equation for the 
heat release can be used in all regions of the cavity with only small error 
introduced. 


: The constant c can be found by conservation of energy as shown 


following. 


Consider the following 3-level molecule: 


E2 


Ey 


Eo 


The quantum efficiency n is defined as the ratio of the energy out divided 
by the total energy available, so for the gain/phase segment under 


consideration. 


(no. of molecules) (E, - E,) AP (A25) 


us (no. of molecules) (E, - Ey) * aH + ap 


| 

4 
assume, for this calculation, that (0,0) is at the corner of the sidewall 
! fi 


Where 
AH = (no. of molecules) (Ey - EQ) ° 
the above expression can be inverted to give 
aH -(=) nee) 
with 
AP =fforcey dx’ dy’ 
and 
wine facar dx’ dy’ 
and the NEP. Then 
P : -(x - x’)/Urtr 
4H = cAz oy fax [ steey Je dx’ 
e e e 


-(x = x’)/Ut 
T(x - x’ )Al(xjy)e ) dx’ (A27) 


—~e 
> 
9 
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Where, recall 


| x > x! 


Since the numerical kinetics routine returns information about the wake 


I(x - x! 
“Yo, x < x’ 
So, 
aay P “(x - x!)/Ut 
AH = cAz fe aise I(x - x’Je 
rp p -x” /Ut 
= cAz | dy [« AI (x3 y dx"e 
—_ «@ 
= chz J» fo AI (x4 y)} “(ce (A28) 
Or, 
AH = c(Az) UTAP 
& 
i cee 
| = (Len) 2 
| : ( 7 Urdz t 
{ 


region itself and not just the heat addition terms, this information must 
be the data used. Thus, for the analytical kinetics model, one must find 


woe ce 


the value for the wake integral: 


x 


4 - -(x!l- x") /Ut 
w(x.y) = fae'ncey = c fexfexrarcer.yre 
eo e eo 


2 cf 20 - x) fawtee - x") AI(x”,y)e 


-(x! = x") /Ur 


~ ~ =(x! = x") /Ut 


= x fewrareety fax’ tex - xf) I(xy’- x%e 


ra -(x = x") /Ut 
2 c fexrarce”» I(x - x”) Jdx’e (A29) 


x 


Ere ence eee lane, hss peer eoemany 2 


So, 
x 
-(x - x") /Ur ssa 
w(x,y) = c | dx"Al(x’,y)UT(1 - e ) (A30) 
° 
So, recalling 
EE Sa (A31} 
n UtAz 
And 
1-G from SIMPGG 
al(x",y) = 2(2 : S) ppp 
ie -(x-x") /Ut 
2 {1 -G\/l-n PPD (x4y)(l-e ) 
woo +H (EeSbs4)/ | 
Az \l G n : (A32) 
Now both numerical and analytical kinetics models return the same array, . 


namely the value of the wake integral throughout the cavity. The effect of 
heat release due to lower level depopulation can be calculated without 
regard to the particular kinetics model chosen. The Fuhs effect is calculated 


in the following manner: 


x (T) 
1) }-w(x(I-1)) 
H(I,J) = a frcx.nex = Sess 3 (A33) 
X(I-1) 


Given this average heat release function, the integral along a 
characteristic can be performed. Note that reflection off the sidewalls 


must be included as can be seen in the following diagram: 


1 » NEP 
N 
| \ 
\ | ‘\ 
\ | \ 
| /\” 
4 
e 
ec 
The contribution at ial due to reflection is therefore found by finding 
the total heat released along the characteristic that reflects at P., 
then adding this to that found along PoP, : 
! | 
The phase shift is found using the Gladstone-Dale law 
n2zil +*Co (A354) i 
The phase change d¢ is : 
2" (CA 
2m (An) (dz — (=e Ao . 
dg = oe fan) (42) = 5 ic. 0% (A35) 


This is then added to that of the unloaded density field to establish the total 
phase change at the gain/phase segment. 
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AIV. SOQ CODE ACCESS 


1. SOURCE CODE 


The following listing represents the source code necessary to update 
the SOQ to include the corrections and modifications described on the 


preceding pages. 


i Toe A? Mowh 


STOENT aaRE 
YINSERT GUL.S3] 
S Tr CIRAQE .NELL) GO TU 850 
‘ WMIR s 1.0 
OELTA 0.” 
NUSTF = 0,n 
ART TE (Rh eABN) [RARE 
BOO FORMAT (/2X e LOHOFO2SOHSEO.4 THAne S%e[ Dee NELTA ANDSe 
C Pe HIStFE ser TA ZFWO AN) wMT 2? SET TU Le Pel Oneseeeeaseuys) { 
ssu CONTINUE 
PINSEHT APO OK, 2D 
wT TEC?) (CUCLZ7) e TZ] 6NUn) 


ay 


REwING 7? 
DELETE ANL e7060G0L.795 
PIELER TE BOL AID eudL eT 


TE C{ICAV.EN.))} THAKES] 
IPCI RABE .FQ.0) FO TO 4D) 
WRITE (29) (CFPFL C12) eT LZ2) sed) ok 
REWIND 29 
NO GAY TeleN03d 
oad Cucl) = CFFLCT) 


Ni) AGT [LaleNPTS 3 
nad K(f) = KKCT) 

OTAMH = ARCllelel) 2-5 

CALL APRTM (DTA He VesVe ede ee ede) 

POwsn, 

REAN (29) (CFFL (TZ) eo L221 eNUH) 0X 

REWIND 29 

110 68m [el eNOs 
6gR POW = Pilw 6 Cult) * COwGI(CUII)) 


Pow = POW 9(XK(2) | AN(L)) 292° (NPI S/NBY) 
WRITE (Ae6H7) POW 
oa? FORMAT(SXy® come PUWER IN FEEDRACK NORMED TO UNITY HY Ge 
f) KELSe7e/) 
SQTPOW = SQRT(KOw) ; 
OU 684 [sl eNO 
5d CFFL(T) = CFFLI(IT) / SUTRPOW 
: ®LELETE GIL A22eb0L 423 
TF (TR ARE LENO) GO TO LON2 
TF CTT TC AMIY, NOT.RESTRII GU To a? 
Cc 22O9Veo CALCULATE FILENVALIE PISGHOOBZS 
TF (POwKkK.GT.9200L) PRwR = 0} 
FIG = SARTIL.21 000. She wWK) 
wr ITE (HeROI FIG 
46 FORMAT (20Ke® ESTIMATED ELGENVALUE = @eE1Se72e/) 
Gy TO 3093 | 
of WRITE (6 0Ra) | 
Ys] FORMAT (20x e9F [AST PASSeIweUT PUwEM NUT UNITY eEIG NUT ESTD4%, 7) ; 
G) TO 1902 | 
| 
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1V02 CUNTINVE 
CALL REGAININCT oN TER) 
1003 TF CICEK.EN.0) GU TU 964 
FUeLETE GIL oR 28 «GOL e434 
PORLETE CAVITY.3 
x 2LIefZL Oe Aare) 


PINGERT CAVITY.141 


TF (LRARE NE .O} Gu TO Lay 


PUELETE CAVITY.1Aa3 


199 WALTF (7) (CuCT™) of Zs] eMut) 


PORLETE CORR] .a? 


CALL DENSY (FLAG eARU oALEN © YLENSUCZ SNAA ONY Ae] e [No WNSYMe THARE) 


PINSERT SON77CY1.19 


IF (IR ORE.GT.O) GO TO Le 


PINSERT CAVITY.217 
6) TO 11 


12 CSCL7) sCMPLA(COS (PHIM) eSIN( PHIM) ) 


PUCLETE NENSY.2 


SUBROUTINE DENSY (FLAG eRAU eXLENe YLENGZSt AR ONPX oNP Yo LF oI NeNS YM, 


xX TAARE) 
@LVSERT DENSY.105 
TE (CTRAWF LEQ. 1) RETURN 
IF (L46.G67.0) GO TO le 
wRITE (60) 3) 


17 FORMAT S/1 9X eSHOee = 6 PFLAG 


4a 15% 
RETURN 
12 CONTINUE 
PINSEST LROPIT 2 
OATA [RARESO/ 
PINSERT GOL .467 
Hata NESTPWsS0./ 
*INSEXT LROPI.L 
CUMMON /KARES/S [RARE 
PINSERT GOL.22 
COMMON /RARFS/S [Rake 
SNELETE LROK1L.3 
KR s[SARE ePLOTS 


SNANAMAN 


INSERT GOL .430 

TNIR(GeI CAV) & IRARE 
®INSERT GOL .437 

ReolOTA (Sel Cav) 
DELETE GOL.439 


De [IN VEMSY%eoH a0, 


2P0S SET TO 0.08//) 


TdaRF 1S FLAG FUR LUALEI)s HAREs UR SEMI@RARE CAVITY 
2 0 FOR LOADED KESUNATOR (OEF AULT VALUE) 
= 1 FOR RARE RESUNATOR 
32 ? FOR SFMI@RAKE RESUNGTUM (UINITY GAINSDENSITY PHOaSeE CHANGF) 


(UNITY GATNed PHASE CHANGF) 


X PROPAGATING PARAMETER o[2/% THARE® *o 13) 


PNELETE GOL .446 


KX RESTHTSININ (Ge [CAV) eZLI(ICAV) sZLOCICAV) se fO[R(SeICAV?) 


PORLE TF GOL.S830 


TF(eNUTSINIT) GU TO 22 


DESTPw 2 0.0 
®INSEAT GDL .2465 

NOH = NPTSenPY 
SUELETE “TQWOR, 26 


70 TF(DESTOW.LEQ L060) GY TO 360 


NOd = NPYONPTS 
NOK2 3 NURe2 
C eee FIND INCIDENT PUWER 
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Sith ots Yai igen iialeslibaelaaiaih labn 


a 


TORN? 


Lord 


AARE 


PrwIN = 0.9 
90 355 [221 eNOR2 02 

355 PPWIN = PPwItt ¢ CUR(IZ)®¢2 «© CuR(I1Z°1) 402 
PrwK = PPwIN@®(X(2) 9X61) )%*29(WwRTS/NPY) /1 000. 
IF (NRE GR.EQ. 1 eOR NREG EQ ee) PPwK SPP WK /WNOWOe2 
TRANS = SQRT(OFS[PWs/PRwH) 

C #@@* SCALE THE HEAM TO THE O&SIREW POWER, 

OO 356 [2 2 1 eNNd2 


356 CUR(TZ) = CUR(IZ) TRANS 
WRITE (604010) DEST PwePrwk o TRANS 
4010 FORMAT (/S5Xe4AHTHE FIELD HAS BEEN SCALED TO DESIRED POWERS 
A AXel 2HDEST PW B9G1204/8X 0] CMPOwK 20Gl2e4/ 
X Bel 2HTRANS BeG1l204/) 


DESI Pw = TRANS 
369 TF CARS CANX) LE cO0N0LUU ceAND ec AWS (ANY) 01 F20.09010U) GO TN 7) 

*DELETE Cl10astG.7 

3 ABC (11 e IMIR o%) eARC(L2eIMI 462) sABC(L3IeIMIR«2) NESIPwW) 

TF (DESTPW.FWe0.0) GO 10 ae 
C @90008 SCALE THE FIELD BACK DOWN. 

9U 368 172] NUS 
356 CU(I7) = CUCIZ) /DESIPW 

ART TF (604000) 

a#ugn FORMAT (/SXe30HTHE FIELO RAS HEEN SCAI FD NOWNS) 
24 CONTINUE 


PTDENT JITTER 
FDELETE GOL.20 
K [CUT eML Ts fOKeIUTHe 
X ICNT2GeICNT25 01 CNT2Z6% 
X [TMeICFK oNCT 
YNELETE GOL.26 
XNSMM (21) eRUV IAN) ePHIA( ZO? eRCiURVE (4) 0OSF (4) eo TLT (4) eI CAVZ (2) 
*UelLeETE GNL.260 
90 173 IZEROZ1lee20 
*0ELETE GOL.314 
OO 177 [7FROz1617 
*INSERT SO077CY1.2165 


Cc 2 23) JITTERS THE dhA4 AN ANGLE ANGJIT 

Cc = 26 DUMMY = LINE 260 [5 TEMPURARILY STORED IN JITTE® TFLOW, 
Cc = 25 UUMMY © LINE 250 [TS TEMPURARILY STORED IN JITTER IFLOW, 
c = 26 OUAMY © LINE 260 [5S TEMPURARILY STORED IN JITTER TFLOW, 


INSERT CORRA./ 
DATA LUTReJITANGeJITUIS 7000 000020/ 
SUELETE GOL .795eS0077CY1.107 
¢ 4 les 177 LAL 197 20/4 el/ 22/7 23/4 2e/ 28/ 26s 
X 015061704180 0190 0200 9210 2 36902350 02600250 e260) o IFLUW 
®UELETE GOL .325+S0Q77CY1 2166 
C 4 les 177 LAL 194 20/ els 22/ 234 24/4 PS 24s 
Kole} 7001400190 0200 021023650230 eee e250 e2hN) TFL IW 
*INSERT GUL .763 


NAMELIST /JITTER/S JITANGeJITOIS 
JITANG = THE ANGLE UF JITTER (IN MICRORADTANS) 


JITOIS 2 THE OISTANCE PRUPAGATED SUCH TraT THE JITTER 
S1GmMa 1S JLITANG®UITTOIS®} .Eo4& 


AOANA 
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: oe : a RMT ase ao aia : _— lets 


c 
*OELE 
C of 
Cc 

230 


231 
1336 


223 


235 
240 
250 


een 


PINSE 


@lmse 


PINSE 


L746 
*UrLE 


TE GUL.327 


SHSSHSHSSSHSSSSSHSSHSHSSHSSHSHSSHSHSHHSHSHHSHSHESHSSHSSSEHSHCHHHOHHHSHSHS 


[UTR 2 [yTQe] 
TE CANOTLINIT) GO TO 231 

READ CINeUITTER?) 

ARC (SefUTRel) 2 JITANG®!] .E-6 

ARC (GeTITQIe2) = YITOIS 

SIGXY SABC (6 eT UTRel) PABC (Oe lUTR ee) 

WRITE (601836) ANC (GelLUTRel) eSIGAY 

FORMAT (45H 998@ SEAM JITTER MUDEL CALLED *¢*%e9TD UF VIAs 
XZ3HTION ANGLE (RANTANS) 2 eGleedeBXe%: STM DEVIATION (SIGAY) 2% 
XG12.5) 

MO 233 IJz1eNOA 

US(IJ) & CUR C2®I Jol) ee2 © CUR(2¢I 4) 202 

NX BR(2) © KCL) 

CALL JITRAGIDXeSIGXY) 

OU 235) [Ja] oNOd 

TJIU & 20T 

CuR(TuTvel) s5ORT(US(IU)) 

CUR(TUTY) = 0.09 

CONTINUE 

CONTINUE 


COMTTANe 

fwNAL -s 1 

G0 TA aay 
ml GOL. 32 

FMITVALENCE CUSCE) eCRIC CLI) o (CURC LY eCuddy) 
NIMENSTON IS (16 3R4) »CUR (327? He? 
aT GlLels 

R eUSeCUw 

REAL JITANGeJITUIS 
RT CYCLED.A? 

NO Le TleleNPTS 

eC{T) = XSAVE CIT) 
TF WAL 

SUBROUTINE QUAL (TR4ASE pe LSAVE oe LPL Fe TITLE pA ce ANS HR ORF eo SI GANG) 


PINSERT QUALI] 


PULLE 
a3 
40 


Yate 


48 
70 


VATA PT73.1461992N6/ 

Th IAL. 1N7 

LF (USTiuAnGe] F295 71070 944 
SIGAY2F OS [Gang 

whl TF (He1836) SIGANY 
FORMAT (45H 800° HEAM JITTER MUUVEL LALLED *eeee5TN DEVI As 
KIGMTTUN ANGLE 2 +6 G2U05 ) 
CALL JT THAG(NKSAVE 6 STORY) 
tiMARBYS 

OO 68 Yale NPTS 

Jl atue) ONPTS 

YO ABM [Be] tlBTS 

T72leul 

TF (Umar. GF .US(I2)9G0 TO 4a 
YUMAXBZUS(T7) 

XPEAK ax (1) 

YRPR OK ar (5) 

TAI Lsee17 

CURT LL LZel) s8Oe Ph OUSCiZ)) 
CUR(TZI7) = 04n 

CUNT TNte 

MER BIIMER/S LON. 


PINSERT MAIN, 300 


i 


SUBROUTINE JLTONG(OXeSTUAY?) 
TeIS SUMROUIT INE MOOLTFIES Tee FAR Feu INTENSITY OISTRIMUTTUN 
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ct Ba 


ANANA 


+ iste a M5 
my ya 


MIOFL THE FRFECTS OF what JETTEWS THE OT TTER TS ASSUMMEN TO 
GAUSSTAN, SINCE THE RhOUL TING INTENSTTY IS THE CONVOLUTION OF THe 
UND ESTURREY PaTENSTTY wil Tre GAUSSTANs THE OPERATION IS PERFORME 
Gy THE FET ON FACH FUNCT [Ow ALUNEs MULTIPLYING [HE KMFSULTS. 

PEM watts Tee LNVERSe FET. UVF eh/26/7H, 

LEVEL PeCiuell 

COMPLE A Cuedl 

COMMON PMELT/SCHCT SING) 0U91 33024) 9K (128) owl NPIS ONPY ODRX DRY 

C MMAR, SCG/SEL V7 993 

DEMENSTON Nw) 

NOTA BT s3,1418 9a, 

SPITS 2 NPTS @ NRTS 

Mean sz 7a 


Malet) 


CI(M) = CMPLX (US(M) 9Ueu) 
10 PPw = PReWweuS(M) 
NNO CL) sNePTS 
NND (2) sNPTS 
NAR = 2°NPTSQ 
NPOZENPTS/2 
CALL FOURT (CleNARetNUe]? 
SIGEXP = 2,9(SIGKY ® PI) ee, 
SINEZ(X(INPTS) @X(1))/e0 © UALS ee 
OF S.5/5STDE 
AO 20 JsalenPfs 
YSWS(( Y=] ) eNF) ee 
TF (JeGT NPV?) YS YS ( (UPNP TQ <1) MDF) 2e0 
IRS (Je) ONRTS 
90 20 yJetenPTS 
RSQ2( (Tal) eQF) aap 
{IF (T.GT.NPU2) xXSUS( (L=NYPTS<1) O0F) ooo 
Kalen 
ON CLURI SCL (K) OF AP (SI GEXP® (ASQ 0 1SQ) ) 
CALL FHURT (CIleNARsNNUe ©)) 
Md 30 KK z) eNPTSU 
WS OKK) CABS ICI (RK)) /NBE TSU 
30 PRPWN = PRPWNeUS (KK) 
PWRFAC = PPW/PPWN 
DU 40 MMe] NPTSY 
40 US(MM) = US (MM) RWwREAC 
WRITE (60100) PWHEAC 
10G FORMAT (//5Xe8THE POWER HAS BEEN SCALEN BY A FACTOR UF sGl2eSe 
X SIN SURROUTING ULTRKG.%//) 
ReETUAN 
ENO 


®UELETE MAIN. 60 


NAMELIST/QLOT/TITLE LOL T eOne ISAV es IPHASE pRAR ORE eSIGANG 


®°NELETE MAIN,230 


210 CALL tar CTPHASE ce ISAVe TUL Te TITLE oRER AS eNB RF eSTGANG) 


*INSERT “MAIN. 22 


NATA STGANG /0.0/ 


APPENDIX B 
SOQ USERS GUIDE UPDATES 
JUNE 1978 TO JANUARY 1979 


INTRODUCTION 


This appendix documents those changes made to the initial SOQ code between 
June 1978 and January 1979. The changes incorporated in the code are those that 
have become generally useful for the physical optics simulation problems which 
have been solved using the SOQ code. The users guide updates are also prepared 
to clarify and correct the initial description of the SOQ code as documented 
and delivered to AFWL on 1 March 1978, in the Preliminary SOQ Users Guide. This 
document supercedes previous written material on SOQ code documentation. The 


organization of the SOQ Users Guide Updates is as follows: 


Section BI New Subroutines 


2. Subroutine CPUTIM 
Subroutine LISTER 


us 


Section BIT Code Changes/Correction 


1. Subroutine ZERN 
BI. NEW SUBROUTINES | 
| 


1. SUBROUTINE ZERN 


Zernike polynomial terms give the SOQ code the ability to model mirrors 
with arbitrary surfaces. This subroutine also provides the determination of 


sensitivity of a given system to the level of these Zernike terms. 


a. Relevant formalism -- The Zernike Polynomials are an orthogonal 
set of polynomials used to describe phase front aberrations. The low order i 
terms of this set correspond to the low order Gauss-Seidel aberrations, such 


as piston, tilt, defocus, astigmatism, coma, and clover. A list of these 


polynomials, Z(k), is given in Table B-1. 


¥ 
¥ 


TABLE B1. ZERNIKE POLYNOMIALS 


1.0 3 | (4r* - 3r*) sin2e Ls 
Rcos@ R*cos4e 


Rsiné R*sinde 


2 


ar? - 1 (10R> - 12R° 


+ 3R) cosé 


R7cos20 (1OR” - 12R° + 3R) sine 


R?sin2e (SR° - 4R°) cos3é 


(3R> - 2R) cose (SR° - 4R°) sin3e 


(3R° - 2R) sine R°cos5@ 
R°cos36 R°sinse 
2 6 4 2 
R"sin3e 20r° - 30rn* + 12R7 - 1 
on” - oR? + 1 7or® - 140r° + seor® - 210R* 
car* ~ 3r2) cos2e 252R19 - 630r® + seor® - 210R4 


+ 30R? - 1 


The phase applied is 


24 
ho => 2nP 2, (R, 9) 
k=1 
= Ao(1,9) 
- = R<l 


° 


Rol (B1) 


If the Zernike radius Ry is specified to be zero, it is a flag to set the 


phase identically equal to zero. | 


b. Fortran formalism -- Subroutine ZERN is called by GDL with IFLOW 
= 24. Namelist ZERNS contains the Zernike radius Ry as well as the coeffi- 
cients of the Zernike polynomials to be applied P(I) [I = (1,24). 


318 


aI >, aptene cba Senne SHY ; ps o » i aa 


Due to excessive use of the FRINGE program, one can also input fringe 
coefficients (PFRNG(1)), corresponding to the 24 Zernike polynomials to be 
applied. The PFRNG coefficients are converted to P coefficients in subroutine 
GDL. 


NAMELIST /ZERNS/ RO, P, PFRNG 


Argument List RO, P 
Commons /MELT/ 
Externals None 


IDENT ZRNIKE computer printouts follow. 


Le iT Peay 


OTe bt she Enb 
tet 
wml Te ETT meats 
A [7b ie ION TINO TE Pom, 
ore Th FLT Tew ee? 
c z Ys Awedt vote Te Aa Few RES Poe iTS OF WAVES 
whim ee THR TT Te a oe 
PP MWSr OT 8h lhe ead 
Lee TOA Fae prtrak geek Ve dere hl ew Temas IVGeNIAIN 
OSE 2T GETTER ee 


TOO HPOESHEHHESSEHEH HSH EEHEEFESH SHES SESE HEH OH HOSE ESEHEHELESEH ESOS OS OS 


C AMOL y Cb hae 
Cem HHROHEHHEHESH HE EEHE HHH SHES HERES SEHEHESHHHHSHEHHEHESHEEELE SEH EHOO EY 
Aor [shee Bz rst w- & | 


th Cowl et ff) wy F> Soe 
Fue { eK 2 6 Ny 
: roy Pad Tepes 

Part PTI) Bing 
‘WY Jud L2te dn 

Dow MED atl} 2 ay 
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SUBROUTINE CPUTIM 
Subroutine CPUTIM has been activated for the CDC computer to print out 
the amount of CPU seconds used by the kinetics package, which is driven by 
Subroutine REGAIN. On the Cyber 176 a system routine 
A = Second (B) 


returns the CP time since start of job, in seconds, to both A and B. 
FORTRAN: 
Argument List: 


IT = 100* time since start of program 
Commons None 


Externals None 


IDENT CPUTIM computer printout follows. 
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3. SUBROUTINE LISTER 


~~ 


Subroutine LISTER was activated so that the output of the resonator. 
design program RESDES or an arbitrary file may be read internally and 
reprinted in the output of the SOQ code. LISTER reads an 80~column file, 
designated as Tape K, and reproduces it in the SOQ-designated system output 
file with pagination defined the same as on Tape K. 

FORTRAN: 
LISTER is called anytime IRSDS, is nonzero in namelist START. 
Argument List: 
K (= IRSDS from START) 
= tape number of the file to be replicated 


Commons: None 


Externals: None 


IDENT LISTER computer printout follows. 
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BII. CODE CHANGES/CORRECTIONS 


The code modifications and corrections included in the code are 
described below by their update file name. The reason for the change, the 
structure, and the listing are included below: 


1.  *ID SOQMAP 


This update provides a cross-reference map to the SO0Q79128 code. The 
first section lists each routine in the order of appearance in the SOQ code 
with its commons and externals. Also given is a list of all routines that 
call it. The second section lists every common block in the SOQ code with 
the subroutines possessing that common block. 
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eis *ID ABCMAP 
Current allocations of the ABC (I, J, K) array are presented here for 
ease in future updating. 
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3. *ID PLTFIX 


Ident PLTFIX modifies the printer-plotting package in the SOQ code. 
This new plot package: 


a. Prints DCALC, IMAX, DCALC FLUX along with the location of the 
center of the beam (DRX, DRY) and the bottom of every iso-inten- 


sity plot 
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Prints a blank for every value of intensity less than 0.01 
*UMAX (UMAX is maximum intensity) and puts a border around 
the outside in column 1 of NPTS and row 1 to NPY 


Allows for selective plotting, based on the new namelist 
parameter KPLOT in namelists PLOT and QLOT. 
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4. *ID ADDPRNT 


$ This section of updates was included to add information on intermediate 
Po printout to CAVITY, STEP, GDL, and TILT: 


a. CAVITY - The incoming and outgoing total flux at 


each gain/phase section is now printed. 


b. STEP - At the beginning of STEP, current values for 
DRX, DRY, RAPTR, NREG, and WNOW are printed, and the 
incoming flux calculated. At the end of STEP, mod- 
ified values of DRX, DRY, NREG, and WNOW are printed 
along with the percent flux lost during the propa- 
gation step. This last parameter (percent flux lost) 
indicates how much of the beam has been propagated 
out of the calculation mesh and, therefore, lost by 
windowing in S-space and K-space (Fourier Transform 
Space). 

c. GDL - At the end of any IFLOW the code now prints 
out total DCALC FLUX, DCALC, and the location and 
magnitude of IMAX. 


d. TILT - Subroutine TILT now prints out the mirror 
radius of curvature necessary to remove the beam 
radius of curvature found by TILT. 
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Miwa S biwe®(X (2) 0K (1) ) SAPPFLOATINOTS/NPY) S19006 
LF (dE Ee Le weit Fie) POWS = PO R/WNOW SO? 
@TISehT STEM 2% 
NOr 2 NRT SPRY 
i) 40) Tae) oN 
tf = 2] 
oop PUOea = Pied © CIO(TT eo), ee? © Cyr (TI) *#2 
PNwA = PNWAXLK (2) OK 1d Pee POELOATONMTOSNPYISTONN 
LF ONME eh ive Lew eo NMF SP e 7) POWA SF POWASWHIWee? 
DELP Ss (BOWASA yA) SOW eI AN, 
LE CORT I ee De Meet IRE Ee D) ANI ec TPRNT otk oth) WRITE (Ae 3000) 
XK AXP AL ce IVYRE AL MRE Te WN We DELP 
DINN FORAATI SOX eo ME CTT ING SURAMITINE STEY. CURRENT PROPAGATION®, 
4 & OnWwaMETFERS: 


X JH 9 OCW z®elslrabe 
: RX /HX,eDNY Beal eeee 
g SRL, ONKE see THe 
UC SME, OwhOW PB eelirieewe 


Be S8X PRP QCEUT FLIER LOST s%@ebl7e4) 
PNELETE STRPLIGI 
16 FURMAT (/aXe} 9H STHE MG TNTENSITY eG eS) 
®INGFQAT STEP.255 
wNOW @ 1.0 
LEC TORNT NE 41) SRATTE (he WAH) NERF AL eo YREAL eNWFE G6 eo wNOw NELO 
OTISE WE STEM ANG 
anOe = 1,0 


GoG 


ONELETE GL eHOT eG IL eAGA 
uU4Max 2 U,0 
KMAK & A(]) 
: YMAX 2 X()) 
DQ 78 Jelenry 
Jl 3 (Ue) ) ONRTS 
OO 7R Psy ese>rs 
1Z = Tel 
AYINT & CIUCTZ) OCONUEICUET/D) 
IF CUMAXeGT.LAYINT) 6 TO 7A 
MAX @ XYEINT 
AMAR 2 X(T) 
YAAX @ K(Y) 
74 PPw =z PRWeXYIN’ 
LE INWE GE QSL ON RES Fie) UMAXBUIMAA/WNO WOO? 
UMAAKK & 11MACs/10049, 
FS WANMAAX B SORT (AMALPODOYMAKOO2) 
oe PELE FE GIN. .AOA 
TF (MSTEP NF ol) Sat TE (079) PRWK eNCALUR oe UMAXK o XMAX 0 YMAL eRADMAX 
®NELETF 2) 870 , 
MUIX BweGleeasArne L2HICALC SoFHe2/Hae lan[ Max SoGlreselUXe 
M2NHLOCBTED AL (CAeY) & Celts LeoSelee's] Powel) s 
XY SOAK eGHAT HAN[US eO127.9) 
MYELFTE Gil 47} 
IF CMSTER CEN. 1) WRITE (Ae TTY} HOAMK eNCALCPetIMAXH eo XMAXK eo YMAX eRalMar 
*HFLETE GOL A723 
MUX Seal eeoe/Ane | AHNCALC seoFReA/HAeLCHIMAX BeGleewel My 
MPDHMLOCETEN AT Cee) = CelsLPerelteelslcodelede 
X SGAXeGMAT BaNTUSe Glee) 
PNP LETTE GOL eS1AeGin o51 4 
{GQN4au = 1 
60 TO 999 


TIL? 
PNFLFTE CYCLE 9.293 
TauRPl = 2, @RAanciiw 
LEC TPS 6G 7) ah TTF (4667) RAC IR e TWORLC 
OOFLETR CYCLES.235 
KL oe SAMMHAGE FOONT CHIRVATUIME 2 RANICIIN BeGlLPA eee te CM/ 

JY%Re® NOTE © Tes CumvaTiiwe CAN AF REMOVED wITH A MTRRORe, 
LVM oO NOT ay WANG Belt eeSe IIMS Ze*maANCUR AS DEFINES) GROVE © / 
10x 6® NEGATIVE VANCUW TS A CONVERGING PHASE FRONT delCH@es 
POne® CAN 6 WEMyEI VET A CONVEX (NEGATIVE RANC) “[YRUR 66s 


ane Ma aA 


x e4)) 


Ca mE ease 


5.  *ID SCLPWK 


Ident SCLPWR modifies the IFLOW = 12 section of GDL to allow for 
scaling of the beam to a specific power TRANS. 


[ORNT SCL PW 


Thu SCL Pak 
HO, 
PLUSFOT FOL aha 

MAG @ 1, 
PELE TE GOL. ade Wie 249] 
© eee TETAS .LE elo) THE FLFELNM TS SCALE! BY SaRTCPRANS) 22 aG 
fC Oe TE CTOANS GT] 60) THE FIELD TS SCALED TO THE PuwER "TRANS 
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| eee 


I] PUL = OIC Che METe 3) 
AMals 3 BC (A cAtLTe'2) 
TRANS = BOLD 
FP CTOANS ot b etet) GO Tu Ass 
Wipwa S U,0 
UN JAA L7s) este 
We PNEw Ss RPNEW © CTL) MCOMIGICUCTZ)) 
PNFWw FS PNE eh CCRC PP MK (CL) OO?) O(NPTSSNMY) SP 000, 
PF UNE (oF ise] aU e IWEG ENG 2) PE WSONE W/WNUWOO? 
GAO LF CTREPS LE Le} ONEw = AMAGOOD 
STUANS SS GET (2.9L 1 7 OME w) 
weT Te (60357) BYC(CLoetEL Te) oAsC(AeMEL Sey) 
ST YSERT GN 2a94 
TF (TRPanS.G°.1) [4n4ta] 


6. *ID TBLUM 


Two errors in subroutine TBLUM are corrected by this ident. The 


following listing is self-explanatory. 


[Meer Tate 


OTT Tet iim 
Tay ae 
SUELETE Pal yum .ar 
TF CAAT SL GF ete ts ARTTE (4 e595) AXTAL 
SORE TE Tal vOM.aar 
CVT 2 (9A CHAS AL FNS OMT S (We OeCPOT) 1) #O( 16/35.) 


7. *ID REMSPH 


Ident REMSPH allows the removal of defocus and/or tilt using a call 
to subroutine QUAL, and to continue with this optimized beam. This optimized 
field can be plotted and written to a local file specified by IWRITE. 


IWRITE. FT.0 sets IW = IWRITE 
IWRITE. LT.O sets in IW = -IWRITE and returns to SOQ 
immediately. 


If desired, the non-optimized field can be read in using ISAV = 1 in 
namelist QLOT. 


PF fit May spry 


PTR YT WY Spe 
SAL 
*T Seah via. 
(PCTS Fete) HD TO Mir 
wide S FE RIG@ py 
To s feamt fe 


332 


TE CT wwRITE LTO) in et ml wet Te 
we LTR Clad) (CrCl ape DT aspe Nie) ek ewX slr eN[ Te Save 
GEwTNi Lo 
wi [TE (Reed) T te TPHASE 

SY FOR 4ATOSS eke ®Cis AAS REC, AWITIEN ON UNT Tel 30% FROM QUAL Se 
X @ a4] Tr [Pear sHeelPss) 

TE (CT aRITECGE 9) 8) To ad 

IPyuT = | 

TE (Ke ET GT et) wRITE (He WN) NETETE 

TE (KMPLOT.LGT. 2) CALL FOLOT(KPLOT) 

LE OTSAVE CRD. i) WR ONDE 7S (CODCI RD EX b eOR) OX gINX CUNY 


AF ate 7 
we Trimed | 
40) CONT Livue H 
“a [8 : 
PRL ETE MAING cP? oo 
AON KePLOTSL OND £ 


HEADS eA QE OT) 


*ID CHGNPT 


Ident CHGNPT increases the flexibility of IFLOW = 6 in two ways: 


a. Reoverlap the beam, letting the code find the 
Original DCALC by setting DIBEAM = 0 


Change the number of points in the beam, by interpolation, 
by specifying NEWNPT and NEWNPY. On a subsequent call to 
START, set NNPTS equal to the value of NEWNPT or NPTS will 
be reset to the previous value of NNPTS. 


ca ae ee OTT Da 


OTe. CHOON T 
AOL 
OUELLETTE (Me MY 
150 [Ler Aas twekaA o 1 
TF Cow To TE) tod TH Tey 
bem RVC TN eih Ct 1) 
ANCL Lah Aen) @ atyedtd 
159 Wher) 2 APC Loe TwKAeh) 
OTSEWT CYCLES A 
NPTS 2 NF anor 
NPY 2 LbhwNRy 
! NOM 8s MAVOere TS 
SPugrof CYCLEQ.43 
Ko NF et Tg ah Vg [AP SM 
ME eo IPT asa SE de Y ARE THE OF STRED) NUMBER OF BOLTS 
SAMPL Tae Tae EFELO, 
fafTOyy 2 Tefo Ci Tee Cat, OF FAULTS TO LMIM & le 
TV PS USE TO RENORMALIZE Te AWE Re SOWATUR FEENRACK FIELD. 
UNSEST Ge AS? 
WENDT 2 
Nery = fy 
IMIWGM 2 } 


Ae ae rem Nene Aen ened Dena ae 


a a el a cra Hi ip SS 


evan 


aged 


PTNSFRT Gol ears 
aC Cmebel) s&s VE yew 
AMC (Pe bol) 3s Ae ahey 
ANC ladeleo}) 2 Tepes 
TF OVD PRAM FIT.) ACC LeLel) S CXCAT ORCI) MEL OAT (NPIS) 
Mie EE TE aN gaa? 
Newel s AUC (Hele tlo.t) 
New ley f= AeCiPeoelelleoedil 
Tels fe @AC(Aslelio wa) 
TF (NG ate TF OGD) ME whi TSNRTS 
LP ONE AMP YF 1) ME NPY SNRY 
ROFL Ss TM SsNE at Tad, 
POELE TE wt ene a 
a AA [ius Ae sh (NPT 
@NELE TRE GIL Ane 
NO AL MYR Leute sy 
SOLE TE GOL HRA a. 
UM 4 MKa] ob wT 
PORLE TE AMVAAST? i 
NE width os wh VRPT Ore Ww py 
AMNOWSQ S 1, 
LP CMIME Geb dele Me elIWE IGF 227) wd ISQswhUwee? 
PwA Sf ROWA/S aN WS : 
Dy Ss ypwhit J yh! pwS ) | 
YY GAY Pes] oti& dye | 
{ 


9.  *ID MISCFX 
This ident corrects minor errors and adds two parameters to namelist 
START. 


° 


a. *D GDL 384,385 This change in format compacts 
this part of the printout to 80 columns for 4 - 
or fewer struts. 


b. *D GDL 884, 885 This change removes the S from 
column 1 so that output can be put on microfiche. 
It also corrects an error in the BARE updates so 
that the CU field is read in at the end of a con- 
verged iteration. 


c.  *D CORR 1.23, 24 This change removes $ from column 
1 in the output. 


d. *D JITTER .83, 86 This change updates the indices 
in the ABC array which were defined originally in 
reverse order. 


e. *D BARE .11 This change corrects the size of the 
loop from MUT to NOB. 


*D Cycle 9. 119, 120 Previously for IPS = 2, 
the iteration counter KOUNT was not updated. 


*I Cycle 9.99 Focal = 1.E50 defaults the radius of 
curvature of the beam to "infinity." 


*I STEP .40 This change activates the IIPS # 0 
option in STEP. Setting DELZ = 0 allows the 


removal of tilt and/or calculated sphere without 


propagating the field. 


*D GDL .827 This statement was redundant. 


*D BARE .86 The parameter RGAIN allows 
the option of not calling REGAIN at the end 


of an iteration. 


The parameter IFLGAP is included so that aperture loads are printed for all 


apertures in the optical train. 


the atl «ESCHER 


*TE ER if “TSCA 


THIS PV Cm E OT AP ee FRG ALA ANY sy MAR AMETERS Tr) 
vaVEL IST STAWT, 


Gf, 
She LE TE to Wwe tg SS 
TNO, IF SPR ITS ee LreF hal At KV CRANE Set Lie Selden] ose AXe 
PV Ag TA AR TR ee TD oe he PTHTHE TASS eAULN SY) 
POPE ETE GOL ieee ZO GRY 
FAD EN AATIS/ Be LEG 44) 77 33N ITERATION [5S CONVEHGED ¢ 
‘ mH NMP friveLaelor [TRAATIANS FAVE e tp LAL) ss) 
ARAN CIP COHN T ZA) oe PZB oititPR) eX epi AersiWY 
WRa TN % 


LISTARO 
HELE TE C RR] AIC HH] 07% 
SEC eAMCAw eT e e LOC T Dobie CTA) oe LOC bers) oe LUC TH) oe LUC LMS) 0 
I SVCD em) OLE La) ob te eS 7 CIM ee OK OM (LMT? JasM/H91)) SK 


ar 


PORPLE TR GT Tew em te JP PPR eee 
ARC CL eLUTo ead = UP TANGE] Fen 
AAC CA] ITH es) =F ITTITS 
PAL STety Sf AMC CLs Luteo) MARE (Ae fl IT em) 
weep hr (mel tay Ad Le LITK em) ea LGry 


Cavity 
SOREL RTE WAawF .] ) 
139 wT Te (7) 


cee qe Te AE 


COHOT2Z) aL 7st ark Or} 


rict 
OUFLETE CYCLE MeL Tee CYCLE 4.120 
eS <a st eto] 
TR CFOS .E L.A) HT) Se 
OL ISka PT CCl bayou 


CCA = Pek Ss 
STF 
PP uSFOL STEP, ou ; 
we ANClin S €2 1C 48 Ba 
CF CDE Zeb tee) RET Sey ‘ 
SOb 
PPELFTE GML eae7 
POELETE RAWF RA 
| IF O98GR INDY CAL SEGATSN(SCT .UITEO) 
MAT 
. OIL VSREST malwel 4 
‘ LOGICAL #GATN 
= PINSFRT AARP OY 
; c T2S9S [5 Tae TARE WhisReR OF TRE AN@CULUMN FILE TO BE COPTED Ta 
) i ceaTw IP Ay LISTE, fF [eSnSsd6 LISTER [5 “OT CALLED, 
c RGOAT™ = oh alLSFe TURNS OFF THE CALL TO wEGaAIN IN TFLOWws7, 
PLUSFYT 4A Na} Se 
Wisqgh s eTwur, 
IPSS = 4 
TFLGAP = 6 
SINSR OT ANINS EAL 
TFL sar = } 
MIELE TE aAdwh 5 
Ro LAAWE PLOTS eRGAINe TRSNS 
@TISFRT AIATNS7 
CMAN SSVTVY 4S RGAT Me TFL GAP 
AL ‘4 
@T sSFoOr Sorel 7 , 


COMMON SOV TYMS HGALNSC TFLGAY 
PELE TE GOL. KH} 
LECT CUTE cE de Le AMM TEL GAM SEI.) GO TH 9G 


10. *ID FXQUAL 


The quality program has been updated to include more options and more 
printouts. See also *ID PROP and *ID RMVSPH for other additions to QUAL 


a. IPRNT. This parameter was added to suppress the 
additional STEP output (from *ID ADDPRNT) when STEP 
! is called from subroutine QUAL. It was also added 
to namelist PROPGT for the same purpose. 


b. The output of the focal plane search was modified to 
: print out more information. 


Ca Additions to QLOT 


RBB (New meaning) 


336 


IWRITE 


PROP 
IRYFF 


KPLOT 


I TABLE 


ICTRD 


(1) 


(see *ID RMVSPHP) 


(see *ID PROP) 


RBB: 


If RBB is input as other than one, 
QUAL will find the quality infor- 
mation for RBB (RA/D) 


IRYTFF.GT.O writes far field 
to unit IRYTFF 


KPLOT.GT.O plots the far field by 
calling IPLOT (KPLOT) 


ITABLE = 0 finds quality table and 
plots information 
= 1 does not do the above 


ICTRD is used for ITABLE = 0, 

= 0 chooses the optimal focal length 
based on the highest 1.0 RA/D quality 
about IMAX, then constructs the quality 
table based on the better of the two 
beam qualities at that focal length 
(Default and same as previous). 

= 1 calculates quality table about cen- 
troid for optimum. 


= 2 calculates quality table about IMAX 
for optimum. 


= 3 finds the optimum value about either 
centroid or IMAX chosen for the highest 
1.0RA/D quality. 


poet F au Abu 


@TORENT FRQUAL 
uae 
@TuSFQT JTTTew. te 
Co ete MRE TS Tee FAW FEED CRNUTERLINE INTENSITY DUE TO A PLANE wave 


c APR ATOHME TO A TAME TER ced wht a CONVERGING LENS UF FOCAL LENGTH 
ot & ABULTED af Tee WE ae FIELD. Toe TT al, PowrR IN THE APEQTUREN 
cr PLANE wAVF ITS Tee Sa4F AS THAT OF THe CURRENS CU FIELD. 


@ENSFOT CYCLE 4.17 
COMMON SSTHOw SZ TRANT 


BML 
CLNSEeT GOLelee 
x 9 | PRNT 

@INSFer GOL,138 
c 
¢ [9eNT TS A FLAG FOW PRINE TYG NREG AND wNOW FROM STFP 
c 20 AONTT OBRINT | 
Cc 2 1 PRINT (OFF AULT) 
@TSEGT wt eS} 

[PHWNT 3 } | 


PINSFAT GOL .S9N 
BBC (Me ISTE Ss 3) 3 TRONT 
INSERT GOL 2593 


{ 
LIPRWE 3 ARC(Se1STER63) © 2901 { 
PINSEQT (ob ela 
COMMON /STPQwI/s TPR 
eek 
@TISFRT JLTTEY. LAA H 
TE CKPLOT SOT. 7) WRITE (6021009) FFTIFL 
WAN FARMAT (LHL een Aass) : 
TE CKOLUT NE.) CALL [PLOT (XPLOT) 
[PenNT z= ] 
IF CTRYTEE Nb ot) wWHeITECIRYTFF) (CUCTK) oe IX2L eNOR) OX eORXODRYONIT e SAVE 
LF CTRYTFF ONE.) APYINU [RYTFF 
TE CTHWYTFE NE.) WRITE (AAD) [RYTFF 
MAA FAWMAT (LOX SPR AR FIELD WAS SEEN WRITTEN TO UNTT@e] a) 
STFR 
STNSERT STEMS 
COMMON /STH wi 7 TPRUT 
Cavity 
eT vSeeTt CaviTy.9 
COMMON sSTRPUwl/ TRPYNT 
#rNSeE WT CAVITY.102 
[Pwiy 2 } 
Qual 
PINSEMT JITTER LII® 
HATA FRTITL 41 4%04 eGHE AM eOHF TEL conf PLesHOTS «2%6H / 
Data NETITL /14%@%4 eOHOPTT GHAI ZF eahi) FL eSHFLI) 02 ¥anm 7 
ATA SAVE 4199067 
[PRNT 2 0 
AAl = 1." 
WHS = 5,7 
@INSFOT Wthdle leo 
C @ee@ [SAVE 2 Qt PENT) [ta FAW FLELO FROM tie [T a, 
@LUSFOT wields 
fc eee [SAVE a2 lo: SAVE NF aw FIELD ON UNIT 7. 


@TNSEIT QUAL 2] 

C #ee [SAVE eel 3s WEAN FAR FLELO FROM UNIT Ye 

@PUSERT CYCLES LAN 

C ee wRITE CH FTEL'Y wETr LENS APPLIEIt (FOC AL LENGTH F) TO UINETT be 
OnFLETE WAL ell? 


CALL PO aww (NETS cd eR eliSe XPE AR 9 YPF AR CRE p PHM) 
TE CISTERP eh deh) CALL POWALW INP TS oD X ek eUS ce XPEAK o YRE AR RNS ePRRS) 
[TF CISTE RP LFO RS CAND CRANE obo?) 
KR CALL PUwWAe (PTS eX eA CUS ce APE AK 0 YPF AR eRHePRRRA) 
PNELFTE QUAL LAL eQiAL.1A7/ 
WRITE (Hel 37) Ot] PK eXCINT eC YCINT eRH] oe PRK cLIMXK gp XPFAK eC YPF AK 
x PWSAVK DR eS TREHL 
TF (RH .NE ele) 
XR ARTE (Ge 137) WH ePRHAWHK eXCINT A YCINT eho ePRRRAK op UMAK o KREAK YE AK © 
x VASAVE NBReSTRERL 
19? FORMAT OS/154 SCALC FLUKA IN eFS 2 0A wL/NSeGlL2e4e® AROUT CENTRO, 
CLAR cL UHCONENTYATE Se ABL2e4/1 SM DCALC FLUX IN F502 06H RL/NZeG1Aebe 
XLOH AAOUT TMad DF ed L2emel 24 CNORNINATES ce 2G612.4/13H TOTAL DCALC « 
XSAPLIIXBOGL? pe Ke 7AM REFERENCE NTAMPTER 2 eF Hee 0/ 
X }94 STHR SL TMTENSETY 20011 04/) 
ANS (WT 2 PRHSK/PwSAVK 100. 
BOSCNT 2 PQASK /PWSAVA e109, 
WRITE (06010) RQOSINT eHUSCNE 
OULD FOMMATISLORS NOTES CENTRINEN anf [Max COORDINATES arf 
© @ CENTIMETEWSe@s/© 8 BADUT [MAX FOR SRL/UBOeGlL2e4e0 
2 10K BY ABMIT TRE CENTROTH) FORM Sxl /N290612.4/) 
PNELETE ‘dtyaL.)2a 
CALL POwwOW (HOTS WK eX eISeKCING oYCINT eRH] oP) 
ZLNSQ = ZENeZIYD 
PINSERT otaLe lo 
IFC TISTEP.WNF em) 50 TY anna 
C eee FIND Bier Ty aR). 
CALL Pl awOW (PTS IX 0K thS oe MPE AK Ce YPF AK eRAS ePRRS) 
CALL PUWWITINPTS eX eh eUS A KCING OVC INT o RAS oe PRAG) 
PHRS = PLNUGE/L1S) 
POWSR 2 PRIAS/ININ. 
PHEAS 3s PWHSOe7L DS) 
PHRSx = PORSsSLYUNN, 
TE (Re OeD 6) GO TO Anu 
FIND BPowFR Tw 2ReRL ZI) 
CALL PiWWOW (NETS eJX ok otiSe XRF AK oe YPE AK eRe PRANTL) 
CALL POWWIAI(NPTS eX eX eS Oo KCINI eo YCINT eRBePRHWA) 
Pele S OVQYNeZI NS) 
PARQHK 2 SORHA/SLONO, 
PQAHYR 2 PUNHWAR71 95.) 
POHPHK 2 PONWA/SINAN, 
C eee WETHION TO CENTIMETERS FOR OUTHUT, 
2900 RCINT = ACTNTOeZLA 
YeIut 2 YCINTe@eZLED 
APE AR = XPF 4K e@7ZLN 
YPGAR 2 YPFAK®7LN 
PH(ISTEP) 2 PAK 
HOTU 2 PRAJPNSAVE*®ING. 
SQCNT 2 PAOR/BNSAVE FIN. 
LF (UROP NE 609) GH TO 340 


POLETE CYCIL& 4 eGae CYCLE 7040 
TE CISTRR EDL) we ETE (Ae 59] U0) 

Sapu FP ywAaAr (S/S 3462 IAM FLuxtet) IN PeRL si ABOUT of 
RC Led THEAL FUCA e AU GHPOTAL emXe ide TAR cENT 
X AXeoe ShwF ty, af 
¥ {ee LEAST HS etree FL ikem Ke TOM (APF AK oe YPF AK) (XCTINTeYCINT) 
NM AMeVHINGE STD Ys 1 he 1G ( ]4e)) ; 
WRITE (05970) [STFU SE ePWSAVK oe PAK CANINT ePHK OHOCNT p STREHL « 
¢ RPE AK eS YRE AM OXCINT A YCINI 

San POMMAT (ISA FellelaweGlAcGeAnehTeAelaelHy eFTeeelHloFaalel4) oAke 

CE ePoi le ae) eA) eth eF2,Gs JIKe lH Cob Me Jet Meek Hoteles 
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Y LKeLA CePA, delete eh Me telh)) 
or FLE TF CYCLE9.17 
Web HM (4) oh 09) PRIA) OXSAVE (12K) eFETITL (20) eNFTITL (20) eSave (10) 
PELE TR CYCLE DAL eCYCLE 3045 
PORTH = =190, 
C @ee FY LUCATIT OM OF MART MEM QUALITY AKROUT [MAX 
NA V79 Lele 
TF (PCL) ott ePrOry) 39 FO 370 
POPPr s R(T) 
sv = 1 
Wa CONTINUE 
C tee FIND LOCATION OF 4a4TMUM QUAL TTY ARDUT TRE CENTROTDe 
Mm) 375 [sales 
TF (PRCT) LE eee re) 3) TO 37S 
POPTH = PRT) 
{SvH = | 
275 CONT Thue 
C eme NETFRMENG FXICAL LENGTH FOW OPTIMAL CALCULATION 
IF CI CTH ch eV eWe (CTRI CE?) TOPTsIS¥v 
TF CICTRE eel) TAMTSISvys 
TF CUCT RI eRe 4d 30) TO JRO 
TORT = TSv 
TF (ROTH GT PORT) [9PT = [SVn 
IAN FORTS kam CT OPT) 
PNELFTE witalel it 
LECT CTR cE Met ee ( LICTIN cE Ve Vee ICT HH cE 0G) eAND ePRALGTSPRA) ) 
x 60 TO 952 
ONFLETE (tal e136 
57 TF (TTantL& Fuel) 30 TO 345 
WRITE (Co eS) ACT NT oO YCINT oF 
Sy FORAAT (S20 0®@ THE FOLLOWING QUALITY ToRALe [5S FOUND AROUTS, 
Ce PYMITMATES (HFA ebvlteeGleese%) FON F Sot 2.4) 
XCINT = KCINTZALEYO 
YCIN?T 2 YCINT/ZLO 
CALL PLTOT INP TS cf eX UMA ee etISe [PL Te 
SORLETE JITTEW.10} 
Sueno IiMe Qral CLPHASE e [SAVE e LPL Te TETL E eRH eC ANS oDReRF oe SI GANG se PROP 
x <PLOTe I wHITEeLTAnL Ee ICTR IRYTFFeNIT) 


MALS 


PNELETE JITTER? 
MAMELIST/QLOT/TI TLE s LOL T eM e LSAV es IPMOSE ohAR ORF CST GANG PROP > 
£ MPLOTe [WRITE eT TAALE ICTR TLHYTEF 

@INGPYT MATN.HO 


c PHIIO s0e PERFORMS FUCAL LENGTH OWTIMEZAT LOM 

c PROP LET CA. CALCINLATRES Qual I Ty FOR TRE NOMINAL FOCAL LENGTH ONLY, 
rc PROP LL Tet. CALCULATES QUALITY FIR TRE CHOSEN FUCAL LENGTH 

c (Fo gs ©PQUP) ONLY. 

c TwH (TR ett. SETS ['¥ 2 TwRITE 

c TwRITE LT. SETS 4a aelywiTF ANt) WETURNS TO SOO IMMEDIATELY. 

C TRYTEF OT. wef TPS THE Fay FIEL) fT) JUNI [RYTFEF 

c RPLOT Gs Ta wMEOTS Tat baw FIELY AY CuLLING [PLOT (KPLOT) 

Cc [TPFaaeF 2 8 Fr WALT TY TARLE a6 PLOTS TNFORMAT TON 

C. 2 1 XOFS NOT DO TRE aROve 

c Icran TS SF) FOR LTARLE = Oe 

¢ = 1 CHQ0SES THE VPTIMAL FOCAL LENGTH 44SEN ON THE HIGHEST 
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11. *ID FILAPR 


Ident FILAPR increases the generality of the aperture routine APRTR by 
adding filleted apertures. 


The outer fillet works by putting a circle in each 


of the four corners of radius Re (input as RFOUT 
or RFMOUT). 


The lightly shaded region is removed by a regular 
rectangular aperture, while the heavily-shaded por- 


tion indicates the region removed by the fillet. 


A central obscuration can be applied in a similar 


fashion. The result is a rectangular aperture with 


rounded corners: 


The input was a list of names for RFIN for 
namelist APTUR and RFMIN for namelist MIROR, 
P Subroutine APRTR has also been modified such that it now prints 


out maximum intensity on both the central obscuration, as well as on 


outer aperture. 


The bare resonator normalization aperture has been generalized to 


include the fillet as well as being any varticular mirror number IM. 
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12. *ID NUDISKT 


Ident NUDISKT modifies the two I/O IFLOWS in GDL, IFLOW = 10 and 
IFLOW = 16. 


a. IFLOW 2 10 Two new options have been added to this 
IFLOW. Multiple fields can now be written to the 
same file by not rewinding the file between writes 
(RWIND = .F.). A file can also be written that can 
read at a terminal (READS = .T.). For this can the 
file is written in the following order: 


TITLE, NPTS, NPY 
(X, [x], I = 1, NPTS) 
DO 141, J 2 1, NPTS 


141 WRITE (IWRITE) ccuft + (J-1) “vers]_ on 
1, NPTS) 


Symmetric fields are unfolded before being written 


to tape for READS = .T. 


b. IFLOW * 16 This IFLOW has been updated so that 
formatted data can be read in as well as written 
out. The format has been modified to include more 


digits. 
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13. *ID CY4KIN 


The capabilities of the numerical SOQ kinetics package have been 
expanded to include oxygen, hydrogen, and R-branch transitions (9.4 y oand). 


a. Oxygen has been upgraded from a structureless mole- 
cule to one that has structure. Therefore, there 
are now kinetics rate equations for the interaction 


of oxygen with the rest of the molecules from the 


b. Hydrogen has been included as a structureless 


{ 
| combustion process. 
} collision partner. 


| 


Cc. Previously the code has used P-Branch transitions 
(10.4 uw). Using GFACT less than 1 now activates the 
9.4 u R-Branch transition. 


In addition to the above major changes, two small additions have been made: 


(1) Input the Gladstone-Dale constant GDC in name- 


list CAVTY2. 
‘4 (2) Account for the gain length by the factor ZFACT, 
* also in CAVTY2. 
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14, *MIRFIX 


The MIRROR subroutine has been modified to calculate the effect of 


power-induced surface curvature when mirror reflectivities other than 
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the design value are used. The parameter, 5, is modified to change the 
center to edge distortion as a function of the mirror reflectivity. The 


parameter, RFLFAC, is used to scale 6 as input through the relation: 


seta (a) 


R = Mirror reflectivity 
P = Incident energy 
d = Design value 


Further, the MIRROR routine has been updated to include the calculation of 
its own value of mirror flux-induced distortion factor when mirrors are 
encountered off axis as noted by PHIAST # 0. This update has not been 
activated, since it would mean input file changes for all users. It is 
included in the code and will be activated by each user, when so desired. 
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15. *ID PROP 


The SOQ code calculation of far field performance is based on the 
analytical equivalence between the Fraunhofer pattern and the propaga- 


tion of a distribution with field curvature, f, a distance Z = f, using the 

Raleigh-Sommerfield formulation of the diffraction integral in the Fresnel 

degree of approximation. The SQQ far field calculation propagates the wave 

distribution, CU, a distance f, determined in a manner which preserves the 

correspondence between near field and far field coordinates, while accurately 
e resolving the energy spectrum in far field coordinates. 


In certain cases, however, it has become necessary to propagate the 
distribution CU to an arbitrary focal plane Z, using the SOQ calculational 
procedure, in order to obtain the effects of beam jitter at a fixed distance 
Z and to obtain the far field information scaled to same focal length. Since 
far field calculations are based on the use of "vacuum" propagation, the far 
field at any plane Z" is simply the scaled distribution at any other plane 
Z'. This can be shown by comparing the far field distributions in terms of 
the Fresnel integrals at two arbitrary focal planes Z' and Z'', where a field 

; curvature of f' = Z' and f" = 2" has been applied to obtain the distribution. ' 
Comparison of these two distributions for the same transmitting aperture size 
leads to the following scaling. ; 


? Zz 
-ik (£ - 2) ik +7 (1 - 3) F 
CU,w = CU, 2 se ze * 5 oy 


i 
f 
! 
: 


And 
Re 2 Re (B4) 


where f is the propagation focal distance obtained in the usual manner from 
the SOQ code, and Z is the "new" scaled propagation distance. 


These changes are incorporated in the SOQ code primarily in subroutine 
QUAL, as documented by the following Fortran changes. 
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16. *ID ECSFIX 


The updates for ECSFIX are included to correct original errors in 
dimensioning Level II variables, and to reduce the resident array sizes at 
load and execution of the code. 
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*ID SEGSOQ 


The SOQ code, as currently configured, is too large to run on the Cyber i 
176 under AFWL Small Core Memory (SCM) restrictions (high speed core). The : 
segmented load option of the CDC NOS/BE loader has been used to reduce : &§ 
execution time SCM requirements without loss of generality of the code. A 
segmentation loader, and the appropriate tree” structure of the code segmen- 


tation is required to take advantage of this feature. 


To incorporate this scheme into the SOQ code, the SOQ code required 
additional 3SLOBAL commons to‘save certain values, as described on the follow- 
ing Fortran listing. A segmentation tree was developed and is listed also. 
Further information on segmentation is available in the CDC/NOS/BE loader 
reference manual. This approach was selected instead of overlay structure 
tecause it is a more powerful tool, even though it is machine specific. 
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18. *WNDOW 


The aerodynamic window subroutine of GDL is used to model the effect of 
an aerodynamic window on the propagated field in a Monte-Carlo sense. The 
aerowindow subroutine simulates a random phase transmission function whose a 
rectangularly distributed random phase information can be selected with 
arbitrary "strength" or variance. This version of AEROW is designed to ae 
simulate the phase field degradation with rectangular probability distribu- 
tion in phase of 0.25). 
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